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Abstract 

In this paper we have studied the flow and heat transfer in a viscous fluid by 
a horizontal sheet. The stretching rate and temperature of the sheet vary with 
time. The governing equations for momentum and thermal energy are reduced 
to ordinary differential equations by means of similarity transformation. These 
equations are solved approximately by means of the Optimal Homotopy Asymp¬ 
totic Method (OHAM) which provides us with a convenient way to control the 
convergence of approximation solutions and adjust convergence rigorous when 
necessary. Some examples are given and the results obtained reveal that the 
proposed method is effective and easy to use. 

Keywords: optimal homotopy asymptotic method (OHAM), film flow, heat 
transfer, unsteady stretching surface. 


1. Introduction 


The flow and heat transfer in a viscous fluid over a stretching surface is a 
relevant problem in many industrial and engineering processes. Examples are 
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manufacture and drawing of plastics and rubber sheets, polymer extrusion, wire 
5 drawing, glass-fiber and paper production, crystal growing, continuous casting, 
and so on. Cooling of stretching surface requires dedicated control of the tem¬ 
perature and consequently knowledge of flow and heat transfer in a viscous 
fluid. Sakiadis [ll, Q, studied the boundary layer flow over a continuous solid 
surface moving with constant speed. Crane analyzed the stretching problem 
10 having in view the fluid flow over a linearly stretching surface. Tsou et al P| 
studied constant surface velocity and temperature. Gupta and Gupta and 
Maneschy et al Q extended the Crane’ work to the stretching problem with 
a constant surface temperature including suction or blowing and to fluids ex¬ 
hibiting a non-Newtonian behavior, respectively. Grubka and Bobba Q studied 
15 the stretching problem for a surface moving with a linear velocity and with 
a variable surface temperature. Wang [sj introduced a similarity transforma¬ 
tion to reduce time-dependent momentum equation to a third-order nonlinear 
differential equation. He analyzed the hydrodynamic behavior of a finite fluid 
body driven by an unsteady stretching surface. The same problem was consid- 
20 ered by Usha and Rukamani Q] for the axisymmetric case. Anderson et al [l(ij | 
analyzed the accompanying heat transfer in the liquid film driven by unsteady 


A5 Q ..a el e, M coneideeea e.ilecee 

and different surface temperature distributions. Vajravelu [13 1 studied the flow 
and heat transfer in a viscous fluid over a planar nonlinear stretching sheet. 
Magyari and Keller [l^ applied the Merkin transformation method to the heat 
transfer problems of stea,dy boundary layer flows induced by stretching surfaces. 
Elbashbeshy and Bazid [l5| studied similarity solution of the laminar boundary 
Iwer equations corresponding to an unsteady stretching surface. Dandapat et al 


16| assumed that the stretching surface is stretched impulsively from rest and 


the effect of inertia of the liquid is considered. The unsteady heat and fluid flow 
has been investigated by Ali and Magyari 171. Liu and Anderson [l^ explored 
the thermal characteristics of a viscous film on an unsteady stretching surface. 
Chen [l^ analyzed the problem of MHD mixed convective flow and heat trans¬ 
fer of an electrically conducting, power-low fluid past a stretching surface in the 
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)resence of heat generation/absortion and thermal radiation. Dandapat et al 


20| studied a thin viscous liquid film flow over a stretching sheet under differ¬ 


ent non-linear stretching velocities in presence of uniform transverse magnetic 
field. Cortell ^ presented momentum and heat transfer for the flow induced 


in a quiscent fluid by a permeable non-linear stretching sheet with a prescribed 
power-low temperature distribution. 

Analytical solutions to nonlinear differential equations play an important 
role in the study of flow and heat transfer of different types fluids, but it is 
difficult to find these solutions in the presence of strong nonlinearity. A few 
approaches have been proposed to find and develop approximate solutions of 
nonlinear differential equations. Perturbation methods have been applied to 
determine approximate solutions to weakly nonlinear problems 22|. But the 


use of perturbation theory in many problems is invalid for parameters beyond 
a certain specified range. Other procedures have been proposed such as the 


Adomian decomposition method 


23|, some linearization methods 2^, 2^ , 


50 various modified Lindstedt-Poincare method s |26|| , variational iteration method 
^ , optimal homotopy perturbation method [ 2 ^^ optimal homotopy asymptotic 
method H - H- 

In the present work we propose an accurate approach to nonlinear differen¬ 
tial equations of the flow and heat transfer in a viscous fluid, using an analyt- 
55 ical technique, namely optimal homotopy asymptotic method. Our procedure, 
which does not imply the presence of a small or large parameter in the equa¬ 
tion or into the boundary/initial conditions, is based on the construction and 
determination of the linear operators and of the auxiliary functions, combined 
with a convenient way to optimally control the convergence of the solution. The 
60 efficiency of the proposed procedure is proves while an accurate solution is ex¬ 
plicitly analytically obtained in an iterative way after only one iteration. The 
validity of this method is demonstrated by comparing the results obtained with 
the numerical solution. 
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2. Equations of motion 


Consider an unsteady, two dimensional flow on a continuous stretching sur¬ 
face, with the governing time-dependent equations for the continuity, momen¬ 
tum and thermal energy B, H, Q, Q: 


du dv 
dx dy 


(1) 


du du dv d^u 


( 2 ) 


dT dT dT d'^T 


( 3 ) 


where u and v are velocity components in the x and y directions, respectively, 
T is the temperature and k is the thermal conductivity of the incompressible 
70 fluid. The appropriate boundary conditions are: 

“0 T ^ rr^ m , "^0 


1 -I- yt ’ 


u — 0, T — Tjx, -l- 


(1 -byt)' 


:( 7 )’ 


at y = 0 


( 4 ) 


U —>■ 0, T —>■ Tn, 


at y —>■ oo 


( 5 ) 


where uq, Tq, Too, 1 are positive constants, c and n are arbitrary and I is a 
reference length. 

If Re = — and Pr = f are the Reynolds number and the Prandl number 
respectively and if we choose a stream function 'I'(a;,y) such that 






u = 


V = -- 


(6) 


dy' dx 

75 then the equation (ED of continuity is satisfied and the mathematical analysis 
of the problems ([2]) and ([3D is simplified by introducing the following similarity 
transformation: 

fiv) 


T ^ 

^ = - 


I \/ife(l + yt)^/^ 


( 7 ) 
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V = 




( 8 ) 




sjy) 

(1 + 7t)'= 


(9) 


Tq being a reference temperature. In this way Eqs. m can be written in the 
form: 



(10) 


(11) 

where prime denotes differentiation with respect to rj. 


Substituting Eqs. (O, ([8|), dU, ([T0|) and (fTTI) into Eqs. 

([2|) and ([3]) , we obtain 

r+//"-r^+A(/'+ir?/") =0 

(12) 

—0" -b /0' - n/'0 -b Afc0 -b ^0') = 0, 
Pr \ 2 / 

(13) 


Here A = ^ is dimensionless measure of the unsteadiness. 

Uo 

85 The dimensional boundary conditions (|3]) and (O become 


^0 ^ 4 - 
u = - - - -Tj (0) at 


I (l+7t)‘ 


Uo 


\/Re{l + 7t)^/^ 


/(O) 


T = roo + 




a;N" 0(0) 


(1 + 7t)^ 


0 = 0 

(14) 

at y = 0 

(15) 

at 2 / = 0 

(16) 


such that for the dimensionless functions / and 0, the boundary/initial condi¬ 
tions become 


fi0) = u, /'(0) = 1, /'(oo) = 0 (17) 
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6 »( 0 ) = 1 , 0{oo) = 0 . (18) 

In addition to the boundary conditions (II3 cLiid. (jlSp . tliG r6C[u.ir6iii6nts 

/(^)>0, %)>0, V77>0 (19) 

must also satisfied [l^ . 


3. Basic ideas of optimal homotopy asymptotic method 

Eqs. (IT^ (or dT^ l with boundary conditions (flTl) (or (IT8l) l can be written 
95 in a more general form: 


N 


($(?7)) 


= 0 


( 20 ) 


where iV is a given nonlinear differential operator depending on the unknown 
function $( 77 ), subject to the initial/boundary conditions: 

d^{r])' 


B 




drj 


= 0 . 


( 21 ) 


It is clear that $( 77 ) = /(? 7 ) or $( 77 ) = 0 ( 77 ). 

Let $o(??) be an initial approximation of $( 77 ) and L an arbitrary linear 
100 operator such as 

d$o(?7)' 


l($o(77))=0, b($o(»?),^^) 


dr] 


= 0 . 


( 22 ) 


We remark that this operator L is not unique. 

If p G [0, 1] denotes an embedding parameter and E is a function, then we 
propose to construct a homotopy [29| - [33|: 


n L(^F{r],p)y H{r],Ci), n(^F{t],p)^ 
with the following two properties: 


(23) 


n 


l(e( 77,0)), iL(7?,a), iv(E(77,0)) 

= l(f(77,0)) =l($o 


( 24 ) 
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n 


L(F(r;,l)), H{r^,C,), A^(F(ry,l)) 


i/(7?,a)iv(ci>(r?)) 


(25) 


where H(r],Ci) ^ 0, is an arbitrary auxiliary convergence-control function de¬ 
pending on variable 77 and on a number of arbitrary parameters Ci, C 2 , Cm 
which ensure the convergence of the approximate solution. 

Let us consider the function F in the form 


F{x,p) = $o(??) +P^’i(? 7 , Ci) C'i) -b ... (26) 


no By substituting Eq. (1261) into equation obtained by means of the homotopy 


n 


l(f(77,p)), H(i^,C,), n(f[i^,p)) 


= 0 


(27) 


and equating the coefficients of like powers of p, we obtain the governing equa¬ 
tion of d>o(a;) given by Eq. (1^ and the governing equation of $ 1 ( 77 , Ci), 
^’ 2 (??,Ci) and so on. If the series (l26l) is convergent at p = 1, one has: 


F(77,1) = $o(?7) + $i(^,C,) + $2(??,C,) + ... (28) 


115 But in particular we consider only the first-order approximate solution 


$( 77 , CO = $o(? 7 ) + $i(??,C 0 , i = 1 , 2 ,..., 777 


(29) 


and the homotopy ((^51) in the form 


n 


l[f{p,p)^, H{ri,Ci), N(^F{r],p)^ 


+p 


i($i(77,CO) -l($o 


H {rj, C, 




(30) 


Equating only the coefficients of p° and p^ into Eq. dani), we obtain the 
governing equation of $o(??) given by Eq. (1221) and the governing equation of 
$ 1 ( 77 , CO i.e. 


B 


L(^$i(77,Ci)) = 77 ( 77 , C0A^(«'o(?7)), 
(^$1(77, Ci), - — -7 = 1,2, ...,777. 


(31) 
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It should be emphasize that ^oiv) and $ 1 ( 77 , Q) are governed by the lin¬ 
ear Eqs. (HH) and m, respectively with boundary conditions that come from 
the original problem, which can be easily solved. The convergence of the ap¬ 
proximate solution (I29l) depends upon the auxiliary convergence-control func¬ 
tion H{r],Ci). There are many possibilities to choose the function H(r],Ci). 
125 Basically, the shape of H(rj,Ci) must follow the terms appearing in the Eq. 
eii). Therefore, we try to choose H{r],Ci) so that in Eq. (IHTl) . the product 
H{r], Ci)N be of the same shape with Now, substituting 

Eq. ((^ into Eq. (l20l) . it results the following residual 


i?(r7,Q) = A^($(77,Q)). 


(32) 


At this moment, the first-order approximate solution given by Eq. (l29ll de- 
130 pends on the parameters Ci, C 2 , Cm and these parameters can be optimally 
identified via various methods, such as the least square method, the Galerkin 
method, the Kantorowich method, the collocation method or by minimizing the 
square residual error: 

J{CuC2,...,Cm)= f R^{v,CuC 2 ,...,Cm) drj (33) 

J a 

where a and b are two values depending on the given problem. The unknown 
135 parameters Ci, C 2 , ..., Cm can be identihed from the conditions: 


dJ 


dJ 


dJ 


= 0 . 


(34) 


With these parameters known (namely convergence-control parameters), the 
first-order approximate solution (1291) is well-determined. 


4. Application of OHAM to flow and heat transfer 

We use the basic ideas of the OHAM by considering Eq. (IT^ with the bound- 
140 ary conditions given by Eq. (II3- We can choose the linear operator in the form: 




(35) 



where > 0 is an unknown parameter at this moment. 

We mention that the linear operator is not unique. Also, we have freedom 
to choose: 




= 


3K 

Kr] + l 




Eq. (|2^ becomes 


^o'"= 0, M0) = fw, $o'(0) = l, $o'(oo)=0. 


which has the following solution 

1 _ p-Krj 

^oiv) = U-\ -^-■ (36) 

145 The nonlinear operator is obtained from Eq. (HU: 

iV/($(ry)) = <l>'"(ry) + iv) - ^'ivf + A($'(r;) + = 0 (37) 

such that substituting Eq. (1551) into Eq. (1571) . we obtain 

= {arj + l3)e~^'^ (38) 

where 

a = ^KA- P = K^-1-KU-A. (39) 

Heaving in view that in Eq. p8p appears an exponential function and that 
the auxiliary function Hf{rj,Ci) must follow the terms appearing in Eq. (l38l) . 
150 then we can choose the function Ci) in the following forms: 

i7/(77, CO = Cl + Car? + (C3 + C4r?)e-^’' + (C5 + ( 40 ) 


or 

i7;(r?, CO = Cl + (Ca + Cgr? + C4r?0e-^’' (41) 

or yet 

CO = Cl + Car? + Car?^ + (C 4 + C^ri)e-^^ + 

+(C 6 + Crr? + C8r?0e-"^'' (42) 
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and so on, where Ci, C 2 , ... are unknown parameters at this moment. 

If we choose only the expression (HOI) for Hf(r], Ci), then by using Eqs. (1551) . 
155 ((5D1) and m, we can obtain the equation in Ci): 


$ 1 '" - = pCi + (aCi + PC2)V + aC2rf 




/3C3 + (oCs + ( 3 Ci)Ti + OLCirj^ 


,-2Kr, . 


/ 3 C 5 + (aCs + PCe)rj ■ 


+aCer]^ 


$i(0) = $i'(0) = $i'(oo) = 0. 


(43) 


The solution of Eq. (|i51) can be found as 


$ 1 ( 77 ) = Ml + 


Ni 


/7aC2 iaCi S/JCa pCi\ 


/ 3cyC2 Q.Ci \ 


1 


+ 


+ 


4it:3 ' 4K^ ' 4K^)^ ' L 108^5 

11/3(74 / llaC^ aCa /3C4\ aC4 2 -2Krj 

~ ^ ~ V 181^4 + ^ J® 


0(7- 




,-Kv . 


850(74 llaCii 


-( 


1150(76 13o(75 1 3/3(76 /3(75 \ 

1728/s:5 288774 288774 24773/® 


(44) 


where 


Ml = - 


3o 


4774 

19o + 1577/3 


277/3 ^ 7a 

Gi- 


4773 


377/3 ^ 5a 
G 2 — 


677/3 


36774 


(7.- 


iVi = 


108775 
3a 


^ 70 + 1277/3 ^ 37a 

G 4 - ^ -^5 ~ 


4277/3 


277/3 7a 

Gi H- 


144774 
377/3 


+ 


4774 
26a + 1277/3 


C 4 


4775 
3a + 477/3 


C 2 + 


864775 
4a + 377/3 


(76 


9774 


C. 


27775 32774 

The first-order approximate solution 
from Eqs. (155)) and (1551) : 


C 5 + 


7a -b 677/3 


Cfi 


(45) 


64775 

for Eqs. m and 013 is obtained 


fiv) = ^iv) = ‘^oiv) + 4>i(r?). 


(46) 


160 III wh&t follows, W6 considGr Eqs. (USD and (fT8)) . In this case, we choose the 
linear operator in the form 


Le[^{v)) =^p" + Kip' 

where the parameter 77 is defined in Eq. ([35]). 


( 47 ) 
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Eq. (1^ becomes 


ipo" + KifQ =Q, 1^0 (0) = 1, </5o(oo)=0. 


Eq. (H5)) has the solution 


165 


Mv) = e-^^. 

The nonlinear operator Ng is obtained from Eq. (USD: 

^e(^‘p{v)'j = ~ + ^(c(p + 

Substituting Eq. (l4^ into Eq. (ICTl . we obtain 

= (mi77 + m2)e“*^’'+ 77136“^^’' 


where 


mi = --KA; m 2 = — - Kf^-1 + cA; m 3 = 1 - n. 

2 Fr 

The auxiliary function Hg(r],Ci) can be choose in the forms: 

Heir], a) =Cr + CsV + (Cg + ^ 1077 ) 6 -^" + (Cn + Ci277)e-"^’' 


(48) 

(49) 

(50) 

(51) 

(52) 

(53) 


or 

H^iv, a) =Cr + CsV + Csr]"^ + {Cw + Cirv)e-^'^ + (54) 

170 or yet 

H;*ir], a) = Cr + (Cs + Cg7])e-^^ + (Cio + Cnv)e-^^^ (55) 


and so on, where C 7 , Cg, ... are unknown parameters. 

If we choose the Eq. (l5^ for Hg, then from Eqs. dSI]), dSSl) and m we 
obtain the equation in Lpi{r], Ci) as 


+ K^pi = m2Cr + {miCr + m2Cs)r] + miCsrf 
m 2 Cg + msCr + (miCg + TO2C10 + msCsir] + miCiorf 
TO3C9 + TO2C11 + (m 3 Cio + rniCii + m2Ci2)r] + miCi2rf 


,-Kn . 


-2Kr) 


-3Kri 


FimsCii + m3Ci2r])e , (j)ii0) = ^i(oo) = 0 . 


(56) 
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Solving Eq. (1551) . we obtain 

2miCs TO 1 C 7 


Viiv) = 


P^ - 


7772^8 \ 2 


V 

TOiCs 3 


iK 
m2CQ 
' 2K^ 


V 


+ 




6-^" + 


msCr 


m2Cs m2C7\ /miCs miCr 
-; —I- In — I : 1- 




K 




2K 


+ 


TmiCio ^ SmiCg ^ 3m2C'io ^ ^ 


( 


4i4:4 

StoiCio 


4^3 4^3 4;;(:3 

miCg 7712(710 , m3Cs\ 


TOi(7io 2 
' 2K^ . 


,-2Kr, , 


2/^2 V 2/^3 2i4:2 2^2 

5(m3(7io + n7i(7ii + TO2(7 i2) 7713(79 + m2(7ii 




6iC2 


+ 


/77l3(7io + 777 i(7ii + r772(7i2 , 5 toi(7i2 \ , 

+l-wP-+ Twf5-j’' + 


771i(7i2 2 


e-3Kv 


flmsCii 777l3(7i2 TO3(7 i2 ^ _4_r-„ j, ^ t \ n 

+ (T2ir3-+ liiKF + lalvDe 7 0i(O)=4i(oo) = O, 


(57) 


175 where 


Pi =- 


^m^Cs 9mi + 6Km2 + 2Km3 


2K^ 4^3 

637711 + 27Km2 + bKnis 


UK^ 


Cg- 


367s:4 


5mi + Pr(TO2 + 37773 ) 

PlO--Pll — 


36P:3 

207772 


77773 


Cl 


(58) 


144P:3 

In this way, the first-order approximate solution (1291) for Eqs. 0 and 0 
becomes 


0{r]) = ip{r]) = ipo{T]) -h (pi{rj,C,). 


(59) 


5. Numerical examples 

In order to prove the accuracy of the obtained results, we will determine the 
180 convergence-control parameters K and Ci which appear in Eqs. 0, 0 by 
means of the least square method. In this way, the convergence-control param¬ 
eters are optimally determined and the first-order approximate solutions known 
for different values of the known parameters fw, A, Pr, n and c. In what follows, 
we illustrate the accuracy of the OHAM comparing previously obtained approx- 
185 imate solutions with the numerical integration results computed by means of 
the shooting method combined with fourth-order Runge-Kutta method using 
Wolfram Mathematica 6.0 software. For some values of the parameters f^, A, 
Pr, 77 and c we will determine the approximate solutions. 
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Example 5.1.a For the first alternative given in the subsection 4.1, we 
consider = —1, A = 1, c = i, n = 1, Pr = 0.7. For Eq. (HSl) . following the 
procedure described above are obtained the convergence-control parameters: 

Cl = -0.0881661632, C 2 = 0.0159074525, C 3 = 101.5499315816, 

Ci = -16.3157319695, = -99.5951678657, Cq = -64.5910051875 

K = 0.7591636981 

and consequently the first-order approximate solution (1461) can be written in the 
form: 


7(?7) = 0.4921333156-t (-0.5668554881-t0.0071536463?7- 
-0.0087518103772 -b 0.00 1 7461596773)e"°-^^®^®^®®®^’’ -b (-0.3980837570 - 
-7.41904137877? -b 2.35914400037;2)e-i-5i832739637, (_o.5271940704-b 

-b6.176450348877 -b 2.3348551362772)e-2-27749i0945,) 

Now, for Eq. (1591) . the convergence-control parameters are: 

C^ = 0.0363993085, Cg = 0.0363993085, Cg = -7.1448075448, 

Cio = 4.3237724702, Cn = 42.1871800319, C 12 = 18.0805214975 

195 and therefore the first-order approximate solution (15911 becomes: 

0 ( 77 ) = (0.3541683003 -b 0.24158769577? - 0.0823906873r?2 -b 
-b0.0060665514773)e“°-^®®^®^®®®^’' -b (-2.6848440528-b 
-b6.366052186677- 1.4238603876772)e"^ ®^®^273963r, 
-b(3.3306757524 - 3.32812400737? - 1.9846972772??2)e“2-27749i0945r) ^g 

In Tables 1 and 2 we present a comparison between the first-order approxi¬ 
mate solutions given by Eqs. dSO]) and (|6T1) respectively, with numerical results 
for some values of variable 7 ? and the corresponding relative errors. 

Example S.l.b In this case, we consider = —1, A = I, c = i, 77 = 1, 
200 Pr = 2. The solution /(??) is given by Eq. (|60|) . The convergence-control 
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Table 1: Comparison between OH AM results given by Eq. Il60t and numerical results for 


fw - 

- 1 , A = 1 




V 

f numeric 

/OHAM> Eq. dlOl) 

relative error = 

/numeric ~ /OHAM 

0 

-1 

-0.9999999999 

1.88 

•10-1® 

1 

-0.1497942276 

-0.1501421749 

3.47 

•10-1 

2 

0.3108643384 

0.3106655367 

1.98 

•lo-i 

3 

0.4604991620 

0.4600705386 

4.28 

•10-1 

4 

0.4887865463 

0.4894195278 

6.32 

•10-1 

5 

0.4919308455 

0.4920394273 

1.08 

•10-1 

6 

0.4921388939 

0.4918417250 

2.97 

•10-1 

7 

0.4921471111 

0.4919782168 

1.68 

•10-1 

8 

0.4921472622 

0.4922151064 

6.78 

•10-® 

9 

0.4921472290 

0.4923440748 

1.96 

•10-1 

10 

0.4921472001 

0.4923640175 

2.16 

•10-1 


parameters for Eq. (l59l) are: 

C 7 = 0.3992391297, Cg = -0.0398233823, Cg = 13.9195482545, 

Cio = -9.8140323466, Cn = -37.1866029355, C 12 = -48.0634410813 

such that the first-order approximate solution (1591) becomes: 

e{r^) = (-2.1207041376- 0 . 056168548877 -h 0.0879369070r;2 - 
-0.006637230377^)6“°-^®®^®^®®®^’' -f e“^-®^®®®’'®®®®’'(7.9716346529- 
-9.1921128520773.231856439677®)e“°-^®®^®®®®®^’' -h (-4.8509305153-h 

-f8.0572360536?7-f 5.275919760677®)e“®-®^^‘‘®i°®’^®’' (62) 

In Table 3 "we present a comparison bet-ween the first-order approximate 
solutions given by Eq. (16211 -with numerical results and corresponding relative 
205 errors. 
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Table 2: Comparison between OHAM results given by Eq. EJ and numerical results for 


fw — 

-1,A = 1, c=i, 

n = l, Pr = 0.7 



V 

n 

^numeric 

^OHAM’ Eq. (Ell) 

relative error = 

(^numeric “ ^OHAmI 

0 

1 

0.9999999999 

8.88 

.10-16 

1 

0.5325816311 

0.5344078544 

1.82 

•lO"® 

2 

0.2137609331 

0.2123010871 

1.45 

•lO"® 

3 

0.0624485224 

0.0627998626 

3.51 

•10-4 

4 

0.0129724736 

0.0141235139 

1.15 

•10-® 

5 

0.0019027817 

0.0018865990 

1.61 

•10-® 

6 

0.0001968219 

-0.0002871987 

4.84 

•10-4 

7 

0.0000144329 

-0.0002516924 

2.66 

•10-4 

8 

8.27 -10-^ 

0.0000468974 

4.60 

•10-® 

9 

1.08 -10-^ 

0.0002281868 

2.28 

•10-4 

10 

7.47 -10-® 

0.0002807824 

2.80 

•10-4 


Example 5.2.a For /„ = 0, A = 1, c = ^ n = 1, Pr = 0.7, the convergence- 
control parameters for Eq. pH]) are: 

Cl = -1.2640611927, C 2 = 0.1680009020, Cg = -34.0575215187, 

Ca = 30.7898356526, C 5 = 37.1281425060, Cg = 13.8590545976, 

K = 1.1203766872 


and therefore, the first-order approximate solution (|46l) can be written in the 
form: 


/(??) = 0.9662722752 -b (1.3563995648 -b 0.0351604322?7 - 
-0.115760505977^ -b 0.0124958644773)e"^-^2°3™®®^2'' -b (-2.5490820161 - 
-1.711111387077- 2.0440805123r72)e"2-^^°^®3^^'‘'‘'' -b (0.2264101761 - 

-0.755240674377- 0.2300192808772)e“3-3®“3°°®^®’^ (63) 
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Table 3: Comparison between OHAM results given by Eq. Il62l l and numerical results for 


fw — 

-1,A = 1, c=i, 

n = 1 , Pr = 2 



V 

n 

^numeric 

^OHAM’ Eq. (ED) 

relative error = 

(^numeric “ ^OHAmI 

0 

1 

0.9999999999 

1.77 

•10-1® 

1 

0.3341908857 

0.3295767993 

4.61 

•10-® 

2 

0.0347540513 

0.0372485103 

2.49 

•10-® 

3 

0.0011305745 

-0.0002320988 

1.36 

•10-® 

4 

0.0000127262 

-0.0002852788 

2.98 

•lO-'i 

5 

-5.44 -10-® 

0.0002991504 

2.99 

•10-1 

6 

-9.15 -10-® 

0.0002884740 

2.88 

•10-1 

7 

-7.96 -10-® 

0.0001372878 

1.37 

•10-1 

8 

-7.05 -10-® 

-0.0000295212 

2.94 

•10-® 

9 

-6.53 -10-® 

-0.0001505702 

1.50 

•10-1 

10 

-5.99 -10-® 

-0.0002044072 

2.04 

•10-1 


210 For Eq. (IMl) . the convergence-control parameters are: 

C 7 = -1.6947892627, Cg = 0.2632100295, Cg = -1.2815579214, 

Cio = 2.6268699384, Cn = 15.8897673738, C 12 = 9.5966071873 

and the hrst-order approximate solution (|59p is: 

%) = (0.2987405005-h 1.1383970916r;- 0.4581387031?72-H 
-l-0.0438683382?73)e"^-^2°^^®®®^^’' -h (-0.1000278783-h 
-hl.48185608457;- 0.58615775577?2)e“2-2'‘°^^33^‘‘4’' -h (0.8012873778- 

-0.59590661657? - 0.7137932043772)e"^-^®“^°°®^®'' (64) 

In Tables 4 and 5 we present a comparison between the first-order approxi¬ 
mate solutions given by Eqs. (1631) and (1641) respectively, with numerical results 
and corresponding relative errors. 
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Table 4: Comparison between OHAM results given by Eq. Il 63 l l and numerical results for 
/^ = 0, A = 1 

_ _ relative error = 

^ /numeric /OHAM’Eq. dMl) 

I /numeric ~ / OHAM I 


0 

-5.50 -10-21 

4.44 -10-1® 

4.44 

.10-16 

1 

0.6894348341 

0.6894914970 

5.66 

•10-® 

2 

0.9167696529 

0.9166682157 

1.01 

•10-4 

3 

0.9608821303 

0.9609858144 

1.03 

•10-4 

4 

0.9659196704 

0.9659030730 

1.65 

•10-® 

5 

0.9662619960 

0.9661631962 

9.87 

•lO-*! 

6 

0.9662759513 

0.9662663279 

9.62 

•10-® 

7 

0.9662762950 

0.9663395358 

6.32 

•10-® 

8 

0.9662763018 

0.9663501433 

7.38 

•lO-*! 

9 

0.9662763032 

0.9663306688 

5.43 

•10-® 

10 

0.9662763043 

0.9663080318 

3.17 

•10-® 


215 Example 5.2.b For /u, = 0, A = 1, c = ^ n = 1, Pr = 2 the firs- 
order approximate solution (1461) is given by Eq. (1631) . The convergence-control 
parameters for Eq. (l5^ are determined as: 

Cr = 0.4652101281, Cs = -0.0724620728, Cg = 14.4924736065, 

Cio = -12.1643720147, Cn = 0.7277740132, C 12 = -56.0450834796 

such that the first-order approximate solution (|59l) may be written as: 

%) = (-1.4030799687-h 0.104261053577-H 0.0880913412?72 - 
-0.0120770121?73)e"^-^^°3^®®®'^^'' -h (3.1476673319- 
-3.808926105777-h 2.7143486992772)e“^-2'‘°’'®^3^'^'‘’’ -h (-0.7445873632-h 

+ 5 . 1973912 OI 877 -h 4.1686190694772)e“^-3®“^°°®^®'' (65) 

In Table 6 we compare between the first-order approximate solutions given 
220 by Eq. ([55]) with numerical results. The corresponding relative errors are also 
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Table 5: Comparison between OHAM results given by Eq. Il64t and numerical results for 


fw — 

0, A = 1 , c = 5 , n 

= 1 , Pr = 0.7 



V 

n 

^numeric 

^OHAM’ Eq. dH 

relative error = 

(^numeric “ ^OHAmI 

0 

1 

1.00 

2.22 

.10-16 

1 

0.4003127445 

0.4006174190 

3.04 

•10-4 

2 

0.1184290061 

0.1183366534 

9.23 

•10-® 

3 

0.0250358122 

0.0254649591 

4.29 

•10-4 

4 

0.0037386526 

0.0032572108 

4.81 

•10-4 

5 

0.0003935287 

-0.0000242886 

4.17 

•10-4 

6 

0.0000291963 

0.0001165652 

8.73 

•10“® 

7 

1.53 -10-° 

0.0003370019 

3.35 

•10-4 

8 

6.39 TO-® 

0.0003255701 

3.25 

•10-4 

9 

8.43 TO-® 

0.0002261157 

2.26 

•10-4 

10 

6.39 TO-® 

0.0001326393 

1.32 

•10-4 


presented. 

Example 5.3.a We consider /„, = 1, A = 1, c = ^ n = 1, Pr = 0.7. The 
convergence-control parameters for Eq. (1461) are given by: 

Cl = 0.6287723857, C 2 = -0.1379103919, C 3 = -64.6127509553, 

C 4 = 52.1862259014, C 5 = 65.7049797786, Cg = 66.9471031457, 

K = 1.6976766716. 

The first-order approximate solution (HSl) one can put as: 

J(rf) = 1.6119245343-h (-0.1095284867- 0.0145745495?7-h 
0.03810895267?2 - 0.0067695650r/^)e-^-^^^^^^^™’^ + (-0.6841717456-H 
4-0.059017118077- 1.5089U67407}^)e-^-^^^^^^^"^^^^ + (0.1817756979- 

-0.627602263477- 0.4839278208772)6“® °®^°^°°^'^®’' (66) 
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Table 6: Comparison between OHAM results given by Eq. Il65l l and numerical results for 


fw — 

0, A = 

1 , c= i, n 

= 1, Pr = 2 



V 

n 

^numeric 

_ relative error = 

^OHAM’ Eq. (EB tt , 

1 ^numeric “ ^OHAmI 

0 

1 


1 

0 


1 

0.1197281855 

0.1187072441 

1.02 

•10-® 

2 

0.0042249398 

0.0041010249 

1.23 

•10-4 

3 

5.09 

•10-5 

-6.04 -10-6 

5.69 

•10-5 

4 

2.41 

•10-" 

0.0001841967 

1.83 

•10-4 

5 

1.67 

•10-® 

0.0000163574 

1.63 

•10-5 

6 

1.42 

•10“® 

-0.0001452859 

1.45 

•10-4 

7 

1.26 

•10“® 

-0.0001791057 

1.79 

•10-4 

8 

1.15 

•10"® 

-0.0001403309 

1.40 

•10-4 

9 

1.09 

•10“® 

-0.0000887807 

8.87 

•10-5 

10 

1.05 

•10"® 

-0.0000493843 

4.93 

•10-5 


225 The convergence-control parameters for Eq. (l5^ , are: 

C 7 = -2.4317156290, Cg = 0.4199773974, Cg = 3.2538989273, 

Cio = 2.8819189890, Cn = 2.6307320394, C 12 = 0.0503403055 

and the hrst-order approximate solution ([5^ becomes: 

e{r]) = (0.0191058146 -f 1.8994242629?? - 0.721678087977^ -f 
-f0.0699962329773)e"^ ®^^®^®®^^®'' -f (0.9928667494-h 
-f0.371287954377- 0.4243916166772)e“^-3®®3®^3'^3^’7 -f (-0.0119725641 - 

-0.1259716711?7- 0.0024710391?72)e“®-°®3°^™^'‘®'' (67) 

In Tables 7 and 8 we present a comparison between the first-order approxi¬ 
mate solutions given by Eqs. (1661) and (1671) respectively, with numerical results 
and corresponding relative errors. 
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Table 7: Comparison between OH AM results given by Eq. Il66l l and numerical results for 
/^ = 1, A = 1 


_ _ relative error = 

^ /numeric /OHAM’Eq. dMl) 

I ./numeric ~ / OHAmI 


0 

1 

1.00 

2.22 

.10-16 

1 

1.5177074192 

1.5176780223 

2.93 

•10"® 

2 

1.6030516967 

1.6030350430 

1.66 

•io-’5 

3 

1.6114161917 

1.6114348208 

1.86 

•10-'^ 

4 

1.6119056438 

1.6119031833 

2.46 

•10-® 

5 

1.6119228465 

1.6119073012 

1.55 

•io-’5 

6 

1.6119232066 

1.6119136285 

9.57 

•10-® 

7 

1.6119232084 

1.6119199331 

3.27 

•10-6 

8 

1.6119232063 

1.6119229505 

2.55 

■10-7 

9 

1.6119232045 

1.6119240509 

8.46 

•10-7 

10 

1.6119232031 

1.6119243982 

1.19 

•10-6 


230 Example 5.3.b For /„, = 1, A = 1, c = ^, n = 1, Pr = 2 the first-order 
approximate solution for f{rj) is given by Eq. (1661) . 

For Eq. (1591) the convergence-control parameters are given by: 

Cr = -0.3023817542, Cg = 0.0519503347, Cg = -27.1653468182, 

Cio = 23.1898068679, Cn = -15.3543205956, C 12 = 31.0226487950. 

The first-order approximate solution (1591) one retrieves as: 

e{r]) = (0.9205703595 - 0.1921602046?? - 0.048718357877^ -b 
-b0.0086583891?73)e"^ ®^^®^®®^^®'' -b (0.2637846962-b 
-b0.942905336677- 3.4149327807?72)e“^-3®®3®^3'‘3^’’ -b (-0.1843550558- 

-2.098658615977- (68) 

In Table 9 we present a comparison between the first-order approximate 
235 solutions given by Eqs. (1681) with numerical results. The corresponding relative 
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Table 8: Comparison between OHAM results given by Eq. dSTj and numerical results for 


fw — 

1, A = 1, c = n 

= 1 , Pr = 0.7 


V 

n 

^numeric 

^OHAM’ Eq. (EZl) 

relative error = 

(^numeric “ ^OHAmI 

0 

1 

1 

0 

1 

0.2625870984 

0.2626175913 

3.04 -lO"® 

2 

0.0500304293 

0.0500306615 

2.32 -10-^ 

3 

0.0067456425 

0.0067634155 

1.77 -lO"® 

4 

0.0006411532 

0.0006125219 

2.86 -lO"® 

5 

0.0000429270 

0.0000457402 

2.81 -10-® 

6 

2.01 -10-6 

2.08 -10-® 

1.88 -lO"® 

7 

5.14 -10-® 

1.35 -10-® 

1.34 -lO"® 

8 

-1.29 -10-® 

6.14 -10-® 

6.16 -10-® 

9 

-1.33 -10-® 

2.24 -10-® 

2.25 -lO"® 

10 

-1.22 -10-® 

7.13 -10-^ 

7.25 40-^ 


errors are presented. 



Fig. 1 Solutions foHAMiv) given by 
(|60|. dMl) and dMl) 
for different values of fu, 

— numerical solution; 

. OHAM solution 


/ 



Fig. 2 Solutions /oHAAiiv) 
obtained from (IMl) . (|63|) and (IMll 
for different values of fw 
— numerical solution; 

. OHAM solution 
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Table 9: Comparison between OHAM results given by Eq. Il68l l and numerical results for 


fw — 

1, A = 1, c = n 

= 1, Pr = 2 



V 

n 

^numeric 

_ relative error = 

^OHAM’ Eq. dMl) tt , 

1 ^numeric “ ^OHAmI 

0 

1 

1.00 

2.22 

.10-16 

1 

0.0288461240 

0.0286378502 

2.08 

•10-4 

2 

0.0002627867 

0.0004342087 

1.71 

•10-4 

3 

1.00 -10-6 

- 1.91 -lO-"^ 

1.91 

•10-4 

4 

1.21 -lO"^ 

- 1.46 -lO-"^ 

1.46 

•10-4 

5 

1.05 -10-^ 

- 3.96 -10-® 

3.97 

•10-® 

6 

9.39 -10-® 

- 4.54 -lO"® 

4.63 

•10-® 

7 

8.68 -10-® 

1.08 -10-® 

9.96 

•10-7 

8 

7.95 -10-® 

8.82 -10-^ 

8.02 

•10-7 

9 

7.44 -10-® 

3.60 -10-^ 

2.85 

•10-7 

10 

6.91 -10-® 

1.18 -10-^ 

4.89 

•10-® 




Fig. 3 Plots of OoHAMiv) given 
by Eqs. (ICTl) . (|M)) and (1571) 
for A = 1, c = y n = 1, Pr = 0.7 
and three values of fw 
— numerical solution; 

. OHAM solution 


Fig. 4 Plots of OoHAMiv) given 
by Eqs. (1571) . (|55|) and (1551) 
for A = 1, c = y n = 1, Pr = 2 
and three values of fw 
— numerical solution; 

. OHAM solution 


In Figs 1 and 2 are plotted the profiles of f{ri) and velocity profile / ( 77 ) 
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respectively for different values of fw It is clear that the solution f{r]) increases 
with an increase of fw and the velocity decrease with an increase of fw The 
condition f (rj) > 0 for 77 > 0 is satisfied. 

245 In Figs. 3-7 are plotted the temperature profiles given for two values of 
Prandl number Pr = 0.7 and Pr = 2 respectively and different values of fw 
From Figs 3 and 4 it is observe that the temperature 9{r]) decreases with an 
increase of the fw for any values of parameter Pr. 


250 From Figs. 5-7 we can conclude that the temperature decrease with of the 
Prandl number and different values of fw 

0 B 




Fig. 5 Plots of OoHAAtiv) 
given by Eqs. (EB) and (IHB 
for A = 1, c = n = I, 

/u, = — I and two values of Pr 
— numerical solution; 

. OHAM solution 


Fig. 6 Plots of 9oHAM{r]) 
given by Eqs. (I64l) and (l 6 ^ 
for A = 1, c = n = I, 
fw = 0 and two values of Pr 
— numerical solution; 
. OHAM solution 
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0 



Fig. 7 Plots of OohamIii) given by Eqs. (|571) and (IS51) for A = 1, c= 
n = I, /u, = 1 and two values of Pr : — numerical solution; 

. OHAM solution 

From Tables 1-9 we can summarize that the results obtained by means of 
255 OHAM are very accurate in comparison with the numerical results. 

6. Conclusions 

In this work, the Optimal Homotopy Asymptotic Method (OHAM) is em¬ 
ployed to propose analytical approximate solutions to the flow and heat transfer 
in a viscous fluid over an unsteady stretching surface. For three values of the suc- 
260 tion/injection parameter f^, the problem admits solutions which are compared 
with numerical solutions computed by means of the shooting method combined 
with Runge-Kutta method and using Wolfram Mathematica 6.0 software. An 
analytical expressions for the heat transfer for two values of the Prandl number 
are obtained. The solution f(r]) increases with an increase of fw and velocity 
265 decrease with an increase of fw The temperature 9{ri) decreases monotonically 
with the Prandl number and with the distance 77 from the stretching surface. 

Our procedure is valid even if the nonlinear equations of the motion do not 
contain any small or large parameters. The proposed approach is mainly based 
on a new construction of the solutions and especially on the involvement of the 
270 convergence-control parameters via the auxiliary functions. These parameters 
lead to an excellent agreement of the solutions with numerical results. This 
technique is very effective, explicit and accurate for nonlinear approximations 
rapidly converging to the exact solution after only one iteration. Also, OHAM 
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provides a simple but rigorous way to control and adjust the convergence of the 
275 solution by means of some convergence-control parameters. Our construction 
of homotopy is different from other approaches especially referring to the linear 
operator L and to the auxiliary convergent-control function Hf and Hg which 
ensure a fast convergence of the solutions. 

It is worth mentioning that the proposed method is straightforward, concise 
280 and can be applied to other nonlinear problems. 
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